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Abstract — In this paper, the effects slip on the peristaltic 
flow of a Carreau fluid through a porous medium in a two 
dimensional channel under the assumptions of low Reynolds 
number and long wavelength is investigated. The flow is 
investigated in a wave frame of reference moving with velocity 
of the wave. The perturbation series in the Weissenberg number 
was used to obtain explicit forms for velocity field, pressure 
gradient per one wavelength. The effects of various pertinent 
parameters on the pressure gradient and pumping 
characteristics are discussed through graphs in detail. 

Index Terms — Carreau fluid; Darcy number; MHD; Slip 
effects 

I. Introduction 

Peristalsis is a series of wave-like muscle contractions 
that moves food to different processing stations in the 
digestive tract. This principle is used in designing the roller 
pumps which are useful in pumping machinery. For instance 
biomechanical pumps are fabricated to save blood or similar 
fluids from any possible contaminations arising out of contact 
with the pumping machinery while pumping the fluid. Since 
many physiological fluids behave like non-Newtonian in 
nature. Many authors have been studied the analysis of the 
mechanisms for peristaltic transport of non-Newtonian fluids 
[1-4]. Carreau fluid model is a four parameter model. The 
peristaltic transport of Carreau fluid by considering different 
flow models were developed by [5-8]. 

Specifically, the non-Newtonian fluids in the presence of 
a magnetic field are very useful in magneto -therapy. The 
controlled application of low intensity and frequency pulsing 
magnetic fields modify the cell and tissue behavior. 
Moreover, the non-invasive radiological test that uses a 
magnetic field (not radiation) to evaluate organs in abdomen 
prior to surgery in the small intestine (but not always). Hence 
magnetically susceptible of chyme can be satisfied from the 
heat generated by magnetic field or the ions contained in the 
chyme. The peristaltic flows of magneto hydrodynamic 
(MHD) fluid have been studied by [9-12]. 

The investigations of blood flow through arteries are of 
considerable importance in various cardiovascular diseases 
particularly arteriosclerosis. In some pathological situations, 
the distribution of fatty cholesterol and artery clogging blood 
clots in the lumen of coronary artery can be considered as 
equivalent to a porous medium. Reference [13] has studied 
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the peristaltic mechanism of a Newtonian fluid through a 
porous medium. Reference [14] has investigated the MHD 
peristaltic flow of a porous medium in an asymmetric channel 
with heat transfer. Reference [15] have studied the Peristaltic 
motion of a Carreau fluid through a porous medium in a 
channel under the effect of a magnetic field. 

Also flows with slip would be use full for problems in 
engineering, for example flows through pipe in which 
chemical reactions occur at the walls, two phase flow in 
porous slider bearings. The initial work on slip boundary 
condition on the peristaltic flow of a Maxwell fluid in a 
channel was discussed by [16]. The effects of slip and non- 
Newtonian parameters on the peristaltic flow of a third grade 
fluid in a circular cylindrical tube were investigated by [17]. 
Effects of slip and induced magnetic field on the peristaltic 
flow of pseudoplastic fluid were analyzed by [18]. Recently, 
[19] have investigated the slip effects on the peristaltic 
transport of a Jeffrey fluid through a porous medium in an 
asymmetric channel under the effect magnetic field. 

In view of these, we studied the effects slip on the 
peristaltic flow of a Carreau fluid through a porous medium in 
a two dimensional channel under the assumptions of low 
Reynolds number and long wavelength. The flow is 
investigated in a wave frame of reference moving with 
velocity of the wave. The perturbation series in the 
Weissenberg number ( We < 1 ) was used to obtain explicit 
forms for velocity field, pressure gradient per one 
wavelength. The effects of various pertinent parameters on 
the pressure gradient and pumping characteristics are 
discussed through graphs in detail. 

II. Mathematical Formulation 

We consider the flow of an incompressible Carreau fluid 
through a porous medium in a two dimensional planar channel 
with flexible walls. It is assumed that the progressive 
sinusoidal waves propagate along the walls of the channel. 
The fluid subjected to a constant transverse magnetic field. 
Induced magnetic field, external electric field, electric field 
due to polarization of charges, heat due to viscous and joule 
dissipation are neglected. The equation of the channel wall is 
given by 

Y = H(X,t) = a + bcos-(X-ct), (1) 

A 

where b , A, c and a are amplitude, wave length, phase speed 
of the wave, mean-half width of the channel respectively, t is 
the time and (X,F) are the Cartesian co-ordinates. Fig. 1 

represents the physical model of the channel. 
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Fig. 1 . The physical model 

We introduce a wave frame of reference (v, y) moving 

with the velocity c in which the motion becomes independent 
of time when the channel length is an integral multiple of the 
wave length and the pressure difference at the ends of the 
channel is a constant. The transformation from the fixed frame 

of reference (X,y)to the wave frame of reference (x, y)is 
given by 

x — X—ct , y = Y, u — U—c , v = V, p(x) = P(X,t) . (2) 

where (w,v) and (U,V) are the velocity components, p 

and P are pressures in the wave and fixed frames of 
reference respectively. 

The constitute equation for a Carreau fluid (given in [20]) 


is 


T = — 


/?— 1 


r 


(3) 


where t is the extra stress tensor, rj^ is the infinite shear rate 
viscosity, rj 0 is the zero shear rate viscosity, T is the time 
constant, n is the dimensionless power-law index and y is 
defined as 
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here ;r is the second invariant of strain-rate tensor. We 
consider in the constitutive equation (3) the case for which 
rj x - 0 and so we can write 


n — 1 


T = ( 1 + ( r ^) 2 ) 2 f 


(5) 


The Carreau model reduces to Newtonian model for n = 1 

(or) r = o. 

The equations governing the flow in the wave frame of 
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where p is the density, k is the permeability of the porous 
medium, o is the electrical conductivity and B 0 is constant 
transverse magnetic field. 

Introducing the non-dimensional variable defined by 
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where Re is the Reynolds number and 8 is the wave number, 
into the equations (6) - (8) (dropping bars), we get 


du dv _ 

— + — = 0 

dx dy 


Re <7 


( du du' 

U h V 


V 


dx dy 


dp 




dx 


- 8 : 


Jf J 

dx dy K 


( 10 ) 

A 


Da 


+ M‘ 


(m + 1) 


J 


Re 8 


3 f dv dv ) dp 


U h V 

K dx dy J 


dy dx dy Da 


( 11 ) 

- — v. (12) 


where r = -2 


XX 


1 + 


f n — l\ 


v 




y 


dx 


T = — 

xy 


T yy = ~ 2S 


1 + 

r n — I s 

We 2 y 2 

du , 2 Si’ ' 


{ 2 J 


^ dx J 


1 + 


( n — O 


V 


We 2 / 2 




dv 

dy 


r = 



( du^\ 

2 

f du e2 

2 

(dv) 

2 

28 


+ 


+ 28 




^dx j 


L dy dx) 





cr k 

and M - ap e H Q I — is the Hartman number and Da - 


% 


a 


is the Darcy number. 

Under lubrication approach, neglecting the terms of 
order 8 and Re, we get 
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here N = JiM 2 + 

V Da 

The corresponding boundary conditions in wave frame of 
reference are given by 
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From Equations (13) and (14) it is seen that p is 
independent of y . So that (13) can be rewritten as 

dp d 
dx dy 

The volume flow rate q in a wave frame of reference is 

given by 

/? 

q - | udy . (18) 

o 

The instantaneous flux Q(X,t) in the laboratory frame is 

h h 

Q(x,t) = ^udy = | (u + \)dy = q + h . (19) 

0 0 
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The time average flux over one period T 


of the 


V c ) 


peristaltic wave is 

_ j / i 
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III. Solution 


Since Eq. (2.17) is non-linear differential equation, it is 
not possible to obtain closed form solution. So, we seek a 
perturbation solution by considering Wiessenberg number 
We as a small parameter. For perturbation solution, we 
expand u,q and p as 


u = a t) + We 2 ii t + CMWe 4 ) 
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dx dx dx 

Substituting these equations in (17) and in boundary 
conditions (15) and (16), we get 

3.1 System of order We° 
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The boundary conditions are 
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3.2 System of order We 2 



dp x d 2 u x ( n — \\ d 


dx dy‘ 


V 




dy 


( du 0 ^ 


v 


dy 


J 


-N 2 u x 


(27) 


The boundary conditions are 


u x + j3 


du x 

+ 

r n — I s 

r du 0 A 

3 

dy 


{ 2 J 

J 



= 0 at y — h 


du x 

dy 


= 0 at y = 0 


(28) 


(29) 


Solving Eq. (24) and then using the boundary condition 
equations (25) and (26), we get 
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where c x — cosh Nh + )3 sinh Nh . 
and the volume flow rate q 0 is given by 
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From Eq. (3 . 1 1 ), we get 
dp 0 _ q (go + h) N 3 
dx sinh Nh - Nhc\ 


(32) 


3.4 Solution for system of order We 2 

Solving Eq. (3.7) using the Eq. (3.10) and the boundary 
conditions (3.8) and (3.9), we get 
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From Eq. (3.14) and Eq. (3.12), we have 
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The dimensionless pressure rise per one wavelength in 
the wave frame is defined as 


Ap = Jo 


dp 

dx 


dx 


(37) 


3.3. Solution for system of order We 
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Da = 0.1, M= 1, p = 0.1 and n = 0.398. 


IV. Results And Discussions 


In order to see the effects of various parameters like 
Wiessenberg number We , power-law index n , slip 
parameter /? , Hartmann number M , Darcy number Da , 


amplitude ratio (p on the axial pressure gradient 


dp 

dx 


we 


plotted Figs. 2-7. From Fig. 2, it is found that the axial 

dp 

pressure gradient — decreases with increasing We . From 

dx 

Fig. 3, it is noticed that the axial pressure gradient 
dp 

— increases with an increase in n . From Fig. 4, it is observed 
dx 


that the axial pressure gradient 


dp 

dx 


decreases with 


increasing p . From Fig. 5, it is found that the axial pressure 
dp 

gradient — increases on increasing M . From Fig. 6, it is 
dx 


noted that the axial pressure gradient — decreases with an 

dx 

increase in Da . From Fig. 7 it is noticed that, the axial 


pressure 


gradient — increases with increasing (j) 
dx 


In order to see the effects of various parameters like 
Wiessenberg number We , power-law index n , slip 
parameter /3 , Hartmann number M , Darcy number Da , 

amplitude ratio (p on the the time-averaged volume flow rate 
Q we plotted Figs. 8-13. From Fig. 8, it is observed that the 


time-averaged volume flow rate Q decreases with increasing 
We in the pumping region (Ap > 0) , while it increases with 
increasing We in both the free pumping (A/? = 0) and 
co-pumping (A p < 0) regions. From Fig. 9, it is found that, 

the time-averaged volume flow rate Q increases with 

increasing n in the pumping region, while it decreases with 
increasing n in both the free pumping and co pumping 
regions. Moreover, it is seen that the pumping is less for 
Carreau fluid than that of Newtonian fluid (n— >l). From 

Fig. 10, it is found that the time-averaged volume flow rate Q 


decreases with increasing /3 in both the pumping and free 
pumping regions, while it increases with increasing /3 in the 

co-pumping region for chosen A p[< 0) . From Fig. 11, it is 

observed that the time-averaged volume flow rate Q 

increases with increasing M in the pumping region, while it 
decreases with increasing M in both the free pumping and co 
pumping regions. Fig. 12, it is found that the time-averaged 

volume flow rate Q decreases with increasing Da in the 

pumping region, while it increases with increasing Da in 
both the free pumping and co pumping regions. From Fig. 13, 

it is noted that, the time-averaged volume flow rate Q 
increases with increasing (j) in both the pumping and free 
pumping regions, while it decreases with increasing (/) in the 
co- pumping region for chosen A p[< 0) . 



Fig. 3. The variation of axial pressure gradient 


dp 

dx 


with n for <p = 0.6 , 


Da = 0.1 , M =1 , p = 0. 1 and We = 0. 1 . 


V. Conclusion 

In this paper, we studied the influence of slip on the 
peristaltic flow of a Carreau fluid through a porous medium in 
a planar channel with the effect of a magnetic field under the 
assumptions of long wavelength and low-Reynolds number 
assumptions. It is observed that, the axial pressure gradient 
and time averaged flux in the pumping region increases with 
increasing n,M and (p , whereas they decreases with 

increasing We , ft and Da . 



Da = 0.1, M = 1 , We = 0.1 and n = 0.398 . 
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Da = 0.1, We = 0.1, >5 = 0.1 and rc = 0.398. 



Fig. 8 The variation of pressure rise A p with time-averaged volume flow 

rate Q for different values of We with ^ = 0.6 , M= 1 , /? = 0.1 , 
Da- 0.1 and n- 0.398. 




x 


Fig. 6. The variation of axial pressure gradient 


— with Da for d> = 0.6 , 
dx 


M = 1, We = 0.1, /? = 0.1and w = 0.398. 


0 

Fig. 9. The variation of pressure rise A p with time-averaged volume flow 

rate Q for different values of n with (j) = 0.6 , M = 1 , fd = 0.1, 

Da = 0.1 and We = 0. 1 . 



Fig. 7. The variation of axial pressure gradient 


— with <j> for We = 0.1, 
dx 


M = 1 , p = 0. 1 and n = 0.398 . 



Fig. 10. The variation of pressure rise A p with time- averaged volume flow 

rate Q for different values of fd with (/) = 0.6 , M = 1 , We- 0.1 
and n = 0.398 . 
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0 

Fig. 1 1 . The variation of pressure rise A p with time-averaged volume flow 

rate Q for different values of M with ^ = 0.6 , We- 0.1, /? = 0. 1 , 
Da = 0.1 and n = 0.398 . 



Fig. 12. The variation of pressure rise A p with time-averaged volume flow 

rate Q for different values of Da with (j) = 0.6 , We = 0.1 , p = 0.1, 

M - 1 and n - 0.398 . 



Fig. 13. The variation of pressure rise A p with time-averaged volume flow 
rate Q for different values of We with (j) = 0.6 , M = 1 , 
P = 0.1, Da = 0.1 and n = 0.4 . 
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